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The  MILES  system  was  designed  anrf  built  as  a  tactical  laser  weaprn 
fire  simulator  (LWFS).  All  of  the  mathematical  analysis  and  experiments 
performed  during  the  Advanced  Development  (AD)  phase  were  to  show  that 
MILES  would  accurately  simulate  tactical  weapon  fire.  The  use  of  MILES 
for  long  range  marksmanship  gunnery  training  will  require  the  kill  zone 
of  the  MILES  transmitter  t>->  he  very  narrow  and  well  defined.  It  is  tie 
purpose  of  this  study  to  construct  mathematical  models  for  long  range 
laser  propagation  which  will  permit  evaluation  of  the  MIIES  equipment 
for  marksmanship  gunnery  training.  The  study  also  includes  an  accurate 
analysis  of  the  Binary  Union  decoder  used  in  the  MILES  receiver.  The 
Binary  Union  decoder  was  designed  to  increase  the  detection  probability 
of  received  MILES  words  under  very  specific  fading  conditions.  The 
currently  available  analysis  makes  several  assumption'  that  are  compatible 
with  the  short  range  problem,  but  for  long  ranges  ^  note  accurate  analysis 
is  required,  for  if  the  short  range  (200  meters)  ar,  ilvsis  were  to  be  ex¬ 
tended  to  the  longer  ranges  an  overly  optimistic  jyslem  performance  results 


A .  Background : 

Iv.  has  long  In  en  recognized  that  an  opt  ical  beam  propagating  only 
several  meters  through  clear-air  turbulence  will  result  in  a  redistri¬ 
bution  of  its  energy,  leading  to  fluctuations  in  the  beam  intensity 
commonly  referred  to  as  "scintillations".  The  statistics  associated 
with  these  scintillations  have  been  the  subject  of  many  theoretical 
and  experimental  investigations  [ 1 1  -  [12].  Much  of  the  theoretical 


work  dealt  with  approximation  schemes  for  the  solution  of  the  wave 


equation  describing  the  propagation  oi  an  optical  beam  through  a 


turbulent  atmosphere.  While  these  approximations  have  led  to  several 
different  theoretical  descriptions  for  the  intensity  fluctuations,  the 
lognormal  model  has  been  more  commonly  accepted  than  the  others.  For 
propagation  paths  up  to  around  200  meters  (or  under  conditions  of  weak 
turbulence)  the  logarithm  of  the  intensity  of  the  distorted  wave  is 
assumed  to  be  normally  distributed;  hence  the  intensity  itself  is  said 
to  have  a  lognormal  distribution.  A  large  number  of  measurements  have 
been  made  over  the  years  which  tend  to  support  the  lognormal  model  in 
the  weak  turbulence  regime  [9]  -  [12). 

For  distances  beyond  200  meters  (or  under  conditions  of  strong 
turbulence)  the  probability  density  function  (pdf)  associated  with  the 
optical  scintillations  is  still  a  matter  of  great  controversy.  Several 
different  models  have  been  proposed  for  this  saturation  regime  where 
multiple  scattering  effects  are  prominent.  These  models  range  from 
perturbed  forms  of  the  lognormal  distribution  [13],  [14],  to  the  K- 
distribution  [6],  [9].  Although  the  K-distribution  provides  reasonably 
accurate  predictions  for  the  statistical  moments  in  part  ol  the  satur¬ 
ation  regime,  it  is  not  theoretically  applicable  in  areas  of  weak 

2  2 

scattering  where  the  mean-square  fluctuation  <1  >/<I>  falls  below  2 

A  new  mathematical  model  is  proposed  here  which  we  believe  is 
applicable  tor  all  propagation  path  lengths  and  conditions  of  turbulence. 
The  model  takes  into  account  both  the  specular  component  of  the  optical 
beam  and  the  diffuse  component(s)  which  begin  to  be  prominent  when 
multiple  scattering  effects  are  important.  We  believe  the  reason  the 
lognormal  model  fails  to  accurately  describe  the  phenomenon  in  the 
saturation  regime  is  that  it  ignores  the  diffuse  component(s)  of  the 
beam.  Owing  to  the  presence  of  both  of  these  components  of  the  beam, 


1 

I 


the  proposed  model  has  characteristics  of  the  lognormal  distribution 
over  short  propagation  paths  whie  approaching  the  negative  exponential 
distribution  as  the  propagation  jiath  tends  to  infinity.  Recent  measure¬ 
ments  made  by  Parry  and  Pusey  [ 9 j ,  as  well  as  our  own  measurements  over 
propagation  paths  from  200  -  3000  meters,  tend  to  support  this  new 


mathematical  model. 


II.  ATMOSPHERIC  MODEL 


Optical  communication  systems  operating  in  a  clear-air  atmosphere 
are  subject  to  many  variations  in  atmospheric  conditions.  To  study  such 
systems  in  general  would  require  a  very  complex  model  of  the  atmosphere.; 
Frequently  it  is  the  case  that  the  particular  phenomenon  of  interest  of 
the  optical  system  may  be  the  result  of  only  certain  atmospheric  varia¬ 
tions  and  other  atmospheric  conditions  can  be  reasonably  ignored.  This 
appears  to  be  the  case  in  studying  the  intensity  scintillations  associated 
with  an  optical  beam  traversing  a  clear-air  atmosphere. 


A .  The  Refractive-Index  Structure  Parameter : 

One  of  the  most  significant  parameters  of  the  atmosphere  for  optical 
propagation  is  the  index  of  refraction.  At  a  point  r  in  space  and 
time  t  the  index  of  refraction  can  be  mathematically  expressed  by 


n(r,t)  =  nQ  +  n^r.t), 


(2. 1; 


where  =  <n(r,l)>  is  the  average  value  over  some  continuously  monitored 
time  period  and  n}  represents  the  random  deviation  of  n  from  its  average 
value . 

For  optical  frequency  signals  the  refractive-index  fluctuations  are 
caused  almost  exclusively  by  atmospheric  temperature  variations.  These 
atmospheric  temperature  fluctuations  originate  primarily  from  the  heating 
of  the  earth's  surface,  but  then  these  variations  are  further  broken  and 
mixed  by  the  wind  which  gives  rise  to  a  continuum  of  different  scale 
sizes.  This  leads  to  rhe  refracture-index  b( ing  proportional  to  the 
temperature  fluctuations 
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Since  the  time  period  associated  with  transmitting  a  laser  beam 
is  of  t  £  order  of  a  few  minutes  or  less,  we  will  be  interested  in 
ref lactive-index  fluctuations  over  similar  time  periods.  Thu*  ,-’£*  can 
assume  the  "frozen-in"  turbulence  model  wherein  the  index  fluctuate  >ns 
ut  a  point  are  generated  by  different  sized  "blobs",  or  turbulent  eddies 
swept  past  the  point  by  the  mean  wind.  With  this  model  we  need  consider 
only  tne  spatial  structure  of  the  turbulence  so  that  (2,;1)  can  be 
rewritten  as 


n(r)  =  nQ  +  n^r).  (2.2) 

The  turbulent  eddies  have  a  lens-like  effect  on  an  optical  beam 

propagating  through  the  atmosphere  [4],  These  eddies  vary  in  size  from 

£  ,  a  scale  of  only  a  few  millimeters  and  called  the  inner  scale  of 
o 

turbulence,  to  an  outer  scale  denoted  by  Lq.  When  the  propagating  beam 
is  near  the  ground  the  dimension  of  Lq  is  roughly  comparable  with  the 
beam  height  above  the  ground.  These  different  sized  eddies  are  respon¬ 
sible  for  both  diffractive  and  refractive  effects  on  the  optical  beam, 
which  results  in  both  constructive  and  destructive  interference  of  the 
beam.  It  is  the  interplay  of  refractive-optical  phase  shift  and 
diffractive-optical  ray  bending  of  the  light  rays  that  account  lor  the 
observed  amplitude  fluctuations,  and  hence,  intensity  fluctuations. 

Assuming  spatial  stationarity  and  homogeneity  of  the  atmosphere, 
we  define  the  autocovariance  function  and  spatial  power  spectrum  of  the 
refractive-index  fluctuations  by  [2] 


Bn(r)  =  <n1(r1)n1(r1  +  r)> 


sin(kr)  dk 


(2.3' 


and 


QO 

<t>  (k)  =  — (  B  (r)r  sin(kr)  dr,  (2.4) 

n  2/k  '0  n 


which  are  three-dimensional  rourier  transforms  of  each  other.  Because  of 
radial  symmetry,  the  integrals  reduce  to  a  single  integral  in  terms  of  the 
radial  variable  r  and  spatial  wave  number  k. 

Tartarski  [l]  suggested  that  the  power  spectrum  had  the  same  form  as 
the  wind  velocity  spectrum,  viz., 

<t>  (k)  =  0.033  C2  k"11/3  exp(-k2/k2  ),  (2.5) 

n  n  r  m 

o 

where  1  =  5 . 92/£  .  The  refractive-index  structure  parameter  is  a 

measure  of  the  intensity  of  the  ref ractive-mdex  fluctuations  (i.e.,  a 
measure  of  the  intensity  of  the  turbulence).  Although  it  is  often 

referred  to  as  the  structure  "constant",  it  typically  ranges  in  values 

-17  -13 

from  10  or  less,  for  conditions  of  weak  turbulence,  up  to  10  or 

more  when  the  turbulence  is  strong.  Over  short  time  intervals  at  a  fixed 

distance,  for  a  constant  height  of  the  propagation  path  above  the  ground, 

it  is  reasonable  to  assume  is  essentially  constant. 

n  1 


b 


Eq.  (2.5)  is  accepted  as  valid  only  for  small  scales  such  that 
k  >  k  .  For  large  scales  such  that  k  <  k  ,  where  k  =  2n/L  ,  the  approx 
imation  is  poor.  A  better  approximation  for  <t>n(k)  in  this  latter  case 
is  given  by  (7) 

<t>  (k)  =  0.033  C2  (k2  +  1/L  )"11/6  expC-k^k2),  (2.6% 

n  n  o  m 

called  the  modified  von  Karman  spectrum. 
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III.  OPTICAL  INTENSITY  FLUCTUATIONS 


In  this  section  we  describe  several  of  the  parameters  of  interest 
in  optical  scintillations  such  as  log-amplitude ,  log- intensity,  norm- 
alizea  variance,  and  log-intensity  variance .  Connecting  relations 
between  these  terms  and  the  lognormal  model  are  discussed  so  as  to 
present  a  background  for  the  math  model  developed  in  Section  V. 

In  the  absence  of  turbulence  the  field  of  an  optical  beam  can  be 
expressed  mathematically  by 

U  (t)  =  Aoeiu)t,  (3.1) 

where  Aq  is  the  (constant)  amplitude  of  the  field,  u)  is  the  angular 
frequency  of  the  wave,  and  i  =  /-"I.  The  incident  intensity  of  the 
undisturbed  field  is  then 


I0  ^  l°o(t)|  2  =  Ao-  (3-2) 

After  propagating  a  certain  distance  into  a  turbulent  medium  the  field 
becomes 

U(t)  =  Ae1UJtel4>  (3.3) 

where  A  and  <t>  represent  the  distorted  amplitude  and  phase,  respectively, 
induced  by  the  medium. 

It  is  customary  to  define  the  log-amplitude  of  the  field  by 

X  =  log(A/AQ)  (3.4) 
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so  that  A  -  A  e*.  Hence,  (3.3)  can  also  be  written 
o  ’ 


U(t)=AeX 

0 


The  intensity  of  the  field  is  now  given  by 


I  =  |  U(t)j  2  =  aVx 


and  we  define  the  log-intensity  of  the  field  by 


log(I/Io)  =  2X> 


(3.5) 


(3.6) 


(3.7) 


which  equals  twice  the  log-amplitude  of  the  field. 

For  both  plane-waves  and  spherical-waves  the  average  intensity  at 
a  point  can  be  taken  as  its  value  in  the  absence  of  turbulence,  i.e., 


<I>  =  I 


Under  this  assumption  it  follows  from  (3.6)  that 


<e2X>  =  1. 


(3.8) 


(3.9) 


The  usual  measure  of  the  strength  of  the  intensity  fluctuations  is  the 
normalized  variance  defined  by 


2  2  2 
2  <1  >  -  <1>  <1  > 


o  = 


-  1. 


<I>2  <I>2 


(3.10) 
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(3.11) 


A Probability  Distribution  of  Log-Intensity : 

Many  theoretical  investigations  have  concentrated  on  the  random 
fluctuations  in  the  log-intensity  or  log-irradiance.  Under  conditions  of 
weak  turbulence  the  log-intensity  appears  to  satisfy  a  normal  probability 
distribution,  and  hence  the  intensity  itself  is  said  to  satisfy  lognormal 
statistics.  To  theoretically  see  why  this  is  so,  we  consider  the  following 


mathematical  model. : 


Let  us  suppose  the  turbulent  medium  to  be  composed  of  a  large  number  of 
independent  slabs,  oriented  perpendicular  to  the  direction  of  propagation  of 
an  optical  beam,  and  such  that  the  thickness  of  each  slab  is  large  compared 
with  the  outer  scale  of  turbulence  (see  Fig.  3-1).  Over  short  distances  the 
scattering  cone  is  narrow  so  most  of  the  energy  of  the  beam  remains  directed 
along  the  axis  of  propagation.  Off-axis  scattering  terms  are  neglected. 


Fig.  1-1.  Multiple  scattering  model  with  weak  scattering 
within  each  layer. 
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The  amplitude  value  of  the  wave  is  initially  Aq.  After  passing 
through  the  first  slab  the  amplitude  is 


A,  =  A  -  m,  A  =  A  (1  -  m, ) , 

1  o  1  o  o  1  ’ 

where  m^  is  a  "small"  random  proportionate.  At  the  second  slab  the 
amplitude  becomes 


A2  =  Al^  ~  m2^  =  Ao(1  "  "  m2^ 


whereas  in  general  after  n  slabs,  we  have 

n 

A  =  A  n  (1  -  m.),  (3.12) 

no.,  j  ’ 

3  =  1  J 


where  the  itu '  s  are  assumed  independent  of  each  other.  Taking  the  natural 
logarithm  of  (3.12)  leads  to 


Ln  A 

n 


00 

Ln  A  +  I  Ln ( 1  -  m . ) , 
o  J 

J  =  1 


(3.13) 


and  if  n  is  sufficiently  large  (n  -*  »)  we  can  invone  the  Central  Limit 

Theorem  to  declare  that  lim  Ln  A  =  A  is  normally  distributed. 

n-*°  n 


Since  we  have  defined 

X  =  Ln(A/Ao)  =  \  Ln(I/Io),  (3.141 
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it  follows  that  x  is  normally  distributed  with  mean 


<X>  =  %<Ln  ( I  / 1 Q )  > 


and  variance 


o2  =  <x2>  -  <x>2. 


The  pdf  for  x  is  therefore 


p(x)  = 


2  2 

exp{-(x  -  <X>)  / 20  }, 

X 


•/2n  ct 


and  by  putting  2\  =  Ln(I/Io),  we  get 


p(I) 


J _ exp{-lLn(I/Io)  -  2<X>]?8ox2} 


2  Jin  o^I 


as  the  distribution  for  I.  If  we  further  introduce 


aLnI  =  <(Ln(I/I0))  >  *  <Ln(I/Io)>2, 


called  the  log-intensity  variance ,  we  find 


o?  ,  ~  4<x2>  "  4<x>2  =  4o  2. 

Lnl  X 


(3.15) 


(3. 1C) 


(3.17) 


(3.18) 


r  3 . 19) 
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That  is,  the  log-intensity  variance  is  four  tiroes  the  log-amplitude 
variance,  Using  (3.19),  we  now  write  (3.18)  in  the  more  familiar  form 
of  the  lognormal  distribution, 


p(l)  =  e*pKL"  1  -  '*2/2a\aT)  ’  (3.20) 


where  v =  Ln  I  +  2 <\>  =  <Ln  I>  (see  Fig.  3-2) 


B.  The  Short  Range  Approximation: 


If  x  is  normally  distributed,  then 


2)^  -  exp(2<x>  +  2c^) 
X 


and  since  <e^>  =  1  (Eq.  (3.9)),  it  follows  that  <x>  = 
Therefore  the  normalized  variance  (3.11)  leads  to 


o  =  <exp(4x)>  -  1 


=  exp(4<x>  +  8o^  ) 


from  which  we  deduce 


2  o 

o  -  exp(4o  *•)  _  i, 


(3.21) 


In  view  of  (3.19),  we  also  note  that 


2  2 
°  =  eXp(°Ln]  )  ~  l' 


(3.22) 
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The  spherical  wave  theory  (1),  [5]  predicts  that  the  log-amplitude 

variance  is  proportional  to  the  index  of  refraction  structure  parameter 

according  to 
n  ■* 


CT2  =  G.  124  C2  k7/6  LU/6,  L»£2/A;  (3.23) 

X  n  ’o’ 

hence  it  follows  that 


0.496  C2  k7/6 
n 


lU/6 


(3.24) 


which  is  accepted  valid  for  all  ranges.  However,  if  we  restrict  the 
log-intensity  variance  such  that  then 


exp (o' 


Lnl 


)  -  1  +  o 


2 

Lnl 


and  thus  (3.22)  leads  to  the  well-known  approximation 


=  0.496  C2 
n 


k7/6  L11/6 


(3.25) 


That  is,  under  conditions  of  weak  turbulence  the  variance  of  the  intensity 

and  the  variance  of  the  log-intensity  are  approximately  equal.  This 

same  relation  is  no  longer  valid  under  conditions  of  strong  turbulence 

2 

where  experimental  evidence  supports  the  point  of  view  that  o  reaches  a 
maximum  value  in  the  saturation  regime  and  then  decreases  steadily 
toward  a  value  of  unity  in  the  limit. 
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IV.  EXPERIMENTAL  MEASUREMENTS 


Recently  collected  data  tend  to  support  once  again  the  lognormal 

model  for  intensity  fluctuations  over  short  ranges  or  conditions  of  weak 

turbulence,  whereas  data  recorded  under  conditions  of  strong  turbulence 

(including  distances  up  to  3000  meters)  illustrate  that  the  mean-square 
2  2 

fluctuations  <1  >/<I>  can  be  as  high  as  6.;  Neither  the  lognormal  moael 
nor  the  newly  developed  K-distribution  (6),  (9j  are  satisfactory  for 
describing  the  observed  phenomenon  very  well  over  all  ranges  of  the 
mean-square  fluctuations. 

A .  Experiment ; 

The  propagation  experiments  which  were  conducted  to  verify  the 
mathematical  modelling  were  carried  out  at  Kennedy  Space  Center's  Space 
Shuttle  Landing  Facility.  The  Shuttle  Landing  Facility  consists  of  a 
concrete  runway  16,000  feet  long  (A. 9  kilometers)  and  300  feet  wide. 

The  area  around  the  runway  was  cleared  of  all  obstacles  and  vegetation 
for  approximately  a  1000  loot  radius.  The  runway  had  weather  monitoring 
stations  at  both  ends  and  midway  down  the  runway.  Also,  satellite 
weather  pictures  were  received  every  15  minutes  of  the  Shuttle  Landing 
Facility  area  along  with  the  entire  east  central  Florida  area.  This 
allowed  the  research  team  to  plan  and  carry  out  the  experiments  under 
well  documented  conditions. 

The  experiments  were  conducted  from  March  17  to  March  28,  1980. 
Tests  were  performed  in  the  afternoon  and  at  night.  The  daytime  temper¬ 
atures  averaged  75  -  80°  F  while  the  nighttime  temperatures  averaged 
66  -  68°  F.  The  humidity  remained  low  both  during  the  daytime  and 
nighttime  experiments.  The  sky  remained  partially  cloudy  during  the 
entire  set  of  experiments  but  the  winds  varied  considerably.  During 
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the  day,  wind  speeds  ranged  from  0-12  mph  (0  -  5.36  m/s),  while  at 
night  the  wind  was  virtually  still.  Wind  direction  was  usually  parallel 
to  the  direction  of  propagation  and  less  than  half  the  time  perpendicu¬ 
lar  to  the  propagation  path. 

The  optical  beam  was  generated  by  a  Spectra  Physics  Model  120  HeNe 
laser  with  a  wavelength  of  632.8  nanometers  and  an  inherent  beam  diver¬ 
gence  of  about  1  milliradian.  The  laser  power  was  15  milliwatts  with 
an  intensity  fluctuation  of  less  than  1%  over  a  period  of  an  hour.  The 
optical  beam  was  used  directly  from  the  laser  without  optics.  A  tripod 
was  used  to  rigidly  mount  the  laser  transmitter  at  1.4  meters  above  the 
runway  surface.  The  runway  was  extremely  flat  which  allowed  the  beam 
to  remain  at  a  constant  height  above  the  surface  for  ranges  up  to  about 
1000  meters.  At  ranges  beyond  1000  meters  the  earth's  natural  curvature 
brought  the  beam  to  within  ^  meter  of  the  runway  surface  at  some  inter¬ 
mediate  point  of  the  path. 

Two  receiving  stations  were  used  to  monitor  the  statistical 
fluctuation  of  the  beam  intensity.  One  station  was  located  at  a  fixed 
distance  of  183  meters  from  the  transmitter  while  the  second  station 
was  placed  at  several  dfferent  distances  down  range  from  both  the 
transmitter  and  the  first  station.  The  short  range  station  was  used 
to  measure  the  first  and  second  statistical  moments  of  the  fluctuating 
laser  beam.  Extensive  testing  at  this  short  range  station  confirmed 
that  the  statistics  were  indeed  lognormal  at  this  range. 

The  short  range  station  consisted  of  a  solid  state  PIN  diode 
detector,  a  high  gain  trans impedance  amplifier,  filter  (set  at  2000  Hz'), 
amplifier,  and  a  special  circuit  designed  to  compute  the  normalized 
second  moment  of  the  fluctuating  signal.  The  aperature  of  the  detector 
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Fig.  4-1.  System  diagram  of  short  range  station. 
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Fig.  4-2 .  System  diagram  of  long  range  station 


Fig.  4-3.  Laser  transmitter  mounted  on  tripode  1.4  meters  above  runway 
surface. 


Fig.  4-4.  Receiver  at  long  range  station. 
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was  \  millimeter  and  the  detector  was  shielded  from  background  light. 
Periodic  checking  of  the  background  light  and  adjustment  of  the  shield 
maintained  a  high  signal- to-background  ratio.  A  system  diagram  of  the 
short  range  station  is  shown  in  Fig.  4-1. 

At  the  long  range  station  a  photomultiplier  detector  was  employed 
in  conjunction  with  appropriate  amplifiers  and  a  Digital  Equipment 
Corporation  MINC-11  computer  which  digitized  the  signal  and  calculated 
the  statistical  moments  (see  Fig.  4-2).  The  aperature  of  the  detector 
was  set  at  50pm  so  that  the  detector  would  function  as  a  point  detector. 

A  high  gain  transimpedance  amplifier  followed  the  photomultiplier  tube 
(PMT),  the  output  of  which  was  filtered  to  reduce  the  shot  noise  generated 
by  the  signal  in  the  PMT.  The  signal  was  then  amplified  prior  to  the  A/D 
converter.  The  A/D  converter  could  accept  signals  bewteen  +5  volts  and 
-5  volts.  The  quantization  level  of  the  A/D  converter  was  2.5  milli¬ 
volts.  So  as  to  obtain  the  lowest  quantization  noise  level  possible, 
the  signal  was  kept  at  the  highest  possible  level  (without  clipping). 

The  sample  rate  by  the  computer  of  the  analog  signal  was  varied 
between  100  and  4000  samples  per  second.  Occasionally  the  background 
light  was  checked  by  turning  off  the  laser  and  then  having  the  computer 
sample  the  background  light  alone.  The  photographs  in  Figures  4-3 
through  4-7  show  the  actual  experimental  equipment  used  during  the  test. 

B.  Results; 

The  analysis  of  the  fluctuating  optical  signal  was  performed  by 
computing  the  first  five  statistical  moments  of  the  measured  signal. 

The  computer  was  programmed  to  calculate  and  display  the  moments  as 
data  was  being  taken.  Normalized  moments  were  calculated  by  dividing 
the  nth  moment  <In>  by  the  average  value  <I>  raised  to  the  qth  power, 


i.e.,  <In>/<I>n.  This  normalization  process  allows  one  to  discern  the 
inherent  statistics  of  the  scattering  process  independent  of  the  power 
of  the  laser  beam  and  allows  for  a  comparison  of  the  statistics  at 
different  ranges.  The  moments  give  information  about  the  shape  of  the 
probability  density  function.  Comparisons  of  the  cumulative  probability 
curves  can  be  deceiving  in  that  different  cumulative  curves  may  appear 
to  be  very  similar  but  yet  the  distributions  from  which  they  came  may 
differ  greatly  out  in  the  "wings"  of  the  distribution.  It  is  these 
"wings"  that  are  of  most  importance  in  long  range  propagation,  and 
since  the  moments  high-light  the  differences  between  the  "wings"  of 
various  distributions ,  we  have  chosen  them  as  a  means  of  discriminating 
between  distributions  rather  than  relying  on  cumulative'  probability 
curves. 

Fig.  4-8  shows  a  plot  of  the  measured  normalized  moments  <In>/<I>n 

2  2 

for  n  =  3,4,5  as  a  function  of  the  normalized  second  moment  <1  >/<I>  . 

It  is  interesting  to;  note  that  the  normalized  second  moment  takes  on 

values  as  high  as  5  and  6.  However,  these  high  values  seem  to  be 

consistent  with  other  recently  measured  data  (9].  For  a  fixed  value 
2  2 

of  <1  >/<I>  the  data  in  the  figure  gives  the  appearance  of  being 
widely  scattered.  Actually  this  is  due  to  the  fact  that  the  moments 
steadily  increase  until  they  reach  a  maximum  value  and  then  begin  to 
slowly  decrease.  As  they  decrease  the  moments  assume  values  lower 
than  corresponding  values  attained  during  the  time  of  increase  for  a 
fixed  value  of  <I^>/<I>^. 

Fig.  4-8  also  shows  graphs  of  the  lognormal  and  K-distributions 

having  the  same  first  and  second  moments  as  the  data.  For  small  values 
2  2 

of  <1  >/<!>  the  data  does  indeed  substantiate  the  lognormal  model  and 
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the  K-distribution  appears  to  fit  the  data  reasonably  well  for  larger 
v*lucs  of  the  normalized  second  moment.  However,  neither  distribution 
provides  the  looping  effect  previously  mentioned  where  the  moments 
reach  a  maximum  value  and  then  decrease  along  different  curves. 

Data  gathered  during  the  nighttime  hours  at  a  distance  of  3000 
meters  is  shown  in  Fig.  4-9.  The  wind  velocity  was  virtually  zero. 

The  wide  scattering  of  the  data  is  attributed  to  the  fact  that  the 
weather  conditions  led  to  a  state  of  nonstationary  statistics. 


V.  A  HEW  THEORETICAL  MODEL  FOR  INTENSITY  FLUCTUATIONS 


Limitations  of  the  lognormal  model  for  describing  random  intensity 

fluctuations  have  been  widely  recognized  in  the  literature.  Even  the 

newly  developed  K-distribution  is  limited  to  only  certain  conditions  of 

turbulence.  Therefore  it  seems  highly  desirable  to  develop  a  universal 

model  which  has  characteristics  of  the  lognormal  model  under  conditions 

2 

of  weak  turbulence  (o^jCCl)  an<*  characteristics  of  the  K-distribution 

2 

for  conditions  of  strong  turbulence 

The  lognormal  model  emerged  under  the  hypothesis  that  virtually  all 
of  the  energy  of  the  optical  beam  is  directed  along  the  axis  of  propa¬ 
gation  and  ignores  the  multiple  scattering  effects  of  the  off-axis  eddies. 
Under  conditions  of  strong  turbulence  the  contributions  from  the  off-axis 
are  more  significant  than  the  line-of  sight  contribution  so  that  ignoring 
these  scattering  effects  can  no  longer  be  valid. 

Following  along  similar  lines  as  in  deWolf's  [4]  physical  model,  we 
believe  it  is  reasonable  to  argue  about  the  existence  of  two  principal 


components  which  contribute  to  the  field  received  at  a  point  downrange 


from  the  transmitter.  One  of  these  components,  called  the  specular 
component ,  'results  from  the  forward  scattering  by  the  large  eddies  along 
the  propagation  axis  as  described  in  the  lognormal  model  in  Section  III-A. 
The  second  component  which  we  refer  to  as  the  diffuse  component (s)  arrives 
at  the  point  receiver  after  multiple  scattering  by  the  off-axis  eddies 
(see  Fig.  5-1).  Thus  the  received  field  is  described  mathematically 


Similar  models  have  previously  been  suggested  also  by  Jao  and  Elbaum  [8] 
and  Fante  [15]. 


SPECULAR  COMPONENT 
ALONE 

(Weak  Turbulence) 


SPECULAR  AND 
DIFFUSE  COMPONENTS 

(Strong  Turbulence) 


Fig.  5-1.  Multiple  scattering  model. 

Based  upon  physical  considerations  it  can  be  argued  that  the 
amplitude  A  of  the  specular  component  essentially  satisfies  the  log¬ 
normal  probability  distribution  and  that  the  phase  angle  0  is  uniformly 
distributed.  These  assumptions  are  consistent  with  deWolf’s  model. 
However,  deWolf  assumes  that  the  amplitude  R  of  the  diffuse  component(s) 
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obeys  a  Rayleigh  distribution.  Surely  this  must  be  the  case  in  the 
limit  where  essentially  all  of  the  energy  of  the  beam  is  contained  within 
the  diffuse  component(s)  and  the  number  of  terms  contributing  to  this 
component  is  large.  In  intermediate  regimes  we  believe  the  diffuse 
component  is  the  result  of  adding  a  small  number  of  terms  so  that  the 
Rayleigh  assumption  is  not  justified.  Therefore  we  choose  to  make  the 
assumption  that  R  satisfies  a  generalized  Rayleigh  distribution  which 
includes  the  Rayleigh  distribution,  among  others,  as  a  special  case. 

As  for  the  specular  component,  the  phase  angle  <|>  is  assumed  uniformly 
distributed. 

The  generalized  Rayleigh  distribution,  also  called  the  m-distri- 
bution,  was  first  introduced  by  Nakagami  [16]  in  connection  with  his 
studies  of  intensity  distributions  for  rapid  fading  in  long  distance 
propagation.  Much  of  his  work  was  done  during  the  1940* s  and  published 
in  Japanese  journals  which  may  account  for  the  fact  that  some  of  his 
work  has  gone  on  virtually  unnoticed.  The  m-distribution  seems  appro¬ 
priate  here  also  in  connection  with  our  problem  because  of  the  close 
smilarities  of  the  phenomena  and  because  it  leads  to  the  Rayleigh 
distribution  as  a  limiting  case. 

A.  Development  of  the  PDF: 

Assuming  the  received  field  within  a  region  of  turbulence  is  given 
by  (5.1),  the  intensity  of  the  field  is  found  to  be 

I  =  |U(t)(2  =  A2  +  R2  +  2AR  cos(<|»  -  9).  (5.2) 


Although  there  are  several  approaches  that  could  be  used  for  calculating 


the  pdf  for  I,  we  do  so  by  utilizing  the  Hankel  transform  (see  Appendix  A) 
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p(I)  *  h  f zC(z)J  (/i  z)dz 


(5.3) 


where 


/  p(A)J  (Az)dA  (  p(R)J  (Rz)dR 

K  o 


(5.4) 


is  a  characteristic  function,  the  function  J  (.)  is  the  Bessel  function 

0 

of  the  first  kind  of  order  zero  ([17],  pp.  355-434)  and  p(A)  and  p(R) 
denote  the  pdfs  associated  with  the  amplitudes  A  and  R,  respectively. 

Under  the  assumption  that  R  satisfies  the  m-distribution,  we  have 
[16]  (see  Fig.  5-2). 


m  2m- 1  , 

2m  R  -mRz/b 

r(m)  bm 


(5.5) 


2  2 
where  b  =  <R  >  -is  the  average  of  the  intensity  R  and  T(.)  is  the  gamma 

function  (see  [17],  pp.  253-294).  Observe  that  for  m  =  1,  (5.5)  reduces 

to  the  well-known  Rayleigh  distribution.  Thus  we  find 


00  00 
j  p(R)JQ(Rz)dR  =  2  mm  f  R2m_1 

•  '  r(m)b,“  o 


e-mR2/b  jo(Rz)dR 


jFj(m;l;-bz^/4m), 


(5.6) 


where  we  have  .used  the  formula  ([18],  p.  716) 


?  2  „  /V  +  U  +  l  ,  „*V 

I  U  ~aX  T  /„  U  r  - )  6 

I  x“e  J^Bx)dx  -  r, 


r(v  +  1) 


-lFl(^2±L;V+1;“fe")  (5,7) 
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Here  ^F^(a;b;c)  is  the  confluent  hypergeometric  function  ([17],  pp. 
503-536).  If  A  is  assumed  to  be  lognormal,  the  first  integral  in  (5.4) 
cannot  be  evaluated  in  closed'  form.  We  therefore  resort  to  the  power 
series  expansion 


J  (x)  -£  (-iLkJ*pl2± 

(kl)2 


(5.8) 


so  that 


(s p(A)J  (Az)dz  *  V*  llil — -  f  A?^  p(A)dA 
L  °  &0  (V\\2  ’o 


,£  (-l)k  (z/2)2k<A2k>  , 


(5.9) 


2k  2 

where  <A  >  denotes  the  kth  moment  of  A  .  Now  (5.3)  becomes 


o»  k  2k 
k-0  2zk  (k!) 


00 

>_  f  z^k^"^F^(m;l;-bz^/4m)J^(*^  z)dz,  (5.10) 


or  owing  to  the  integral  formula  ([18],  p.  873) 


iLm~  G?i  / 


[x2p. F.(a;b;-Ax2)J  (xy)dz  •  - - 

4  1  1  V  r(a)  y  p+i 


2^  / Y_  ^  \  (5  u) 

23  (  4X  v+1  ,  1— v  ,  J* 

\  —  +  p,  a,  -j~+  P/ 


(5.10)  reduces  to 


°*  /  ,sk  ,2k 

P<X)  -E  —  2  ■ - 
k»0  (k!)  r(m) 


r-(k+l)  21  /  ml  1,1 


23  l  b  k+.l 


i.  1  ) 

1,  m,  k+1/’ 


(5.12) 
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which  we  shall  henceforth  refer  to  as  the  G-dlstributlon  or  universal 

models  The  function  Gnn  (.)  is  a  generalized  function  known  as  the 
— —  P9 

Meljer  G-function  ([19],  p.  206). 

Although  (5.12)  is  a  representation  for  the  pdf  associated  with 
Intensity  fluctuations,  it  does  not  lend  itself  in  its  present  form  to 
easy  computations  because  of  the  complexity  of  the  G-functioh.  In 
Appendix  B  we  relate  this  particular  G-function  to  the  associated 
Laguerre  functions,  but  even  so  the  resulting  expression  is  still  quite 
complex.  Further  refinement  of  this  pdf  is  called  for  in  order  to  make 
it  more  tractable. 


B.  Theoretical  Moments: 

For  purposes  of  comparing  a  probabilistic  theoretical  model  with 
experimental  data  we  believe  a  comparison  of  the  predicted  (normalized) 
moments  with  the  actual  (normalized)  moments  provides  a  better  criterion 
for  validating  the  model  than  does  a  comparison  of  the  distributions 
themselves  or  their  cumulative  probabilities.  With  few  exceptions, 
the  moments  uniquely  determine  the  distribution. 

While  the  pdf  (5.12)  does  not  readily  lend  itself  to  direct  com¬ 
putations  for  the  pdf  or  the  cumulative  probability,  it  does  lead  to  a 
closed  form  expression  for  the  moments,  rnus  if  ye  use  the  integral 
formula  ([18],  p.  897) 


/  lai . V\ 

Hh . b/* 


m  n 

ri  r  (b,  +  p)  II  r  (l  -  a,  -  P) 
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.  n  r  (i  -  b;,  -  p>  m  r(a<  +  p) 

J-mfl  3  3 
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(5.14) 


where  the  series  is  now  finite.  That  is,  all  terms  such  that  k  ^  n  +  1 
reduce  to  r.ero  because  of  the;  factorial  (n  -  k)!  in  the  denominator. 
Also, 


r(k  -  n)  ,  ,%k  ,  ,  ,  .,n  (~l)k  n! 

— ■ a -  *  (-1)  n(n  -l)...(n  -  k  +1)  “  - — i - 


r(-n) 


(n  -  k) ! 


and  thus  (5.14)  is  equivalent  to 


I2  bn~kr(m  +  n-  k)<A2k>  , 

n-k  . 
m  T(m) 


(5.15) 


where 


k!  (n  -  k) ! 


In  order  to  make  comparisons  of  the  moments  of  various  distri¬ 
butions  they  should  be  normalized  so  that  the  mean  value  is  unity. 
To  do  so,  we  first  set  n  =  1  in  (5.15)  to  get 


<I>  =  b  +  <A>*»b(l  +  r) 


(5.16) 
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where  r  »  <A  >/b  *  <A  >/<R  >  is  the  power  ratio  of  mean  intensities  of 
the  specular  to  diffuse  components.  It  is  also  convenient  to  introduce 
the  normalized  moments  [16],  [20] 


r(m  +  k) 
m^  r(m) 


(5.17) 


<A2k>  k(k-l)/2 

ak  *  IT  "  x 

<A  > 


4  2  2 

where  X  ■  <A  >/<A  >  .  Hence  we  now  write 


<In>-bnE  W  Vn-k  ak  r 
k*0  W 


so  upon  division  by  <I>  ,  we  finally  obtain; 


—  „  • - I - £  (k)2  Vk  ^  rK. 

<I>"  (1  +  r)n  t£b'  ' 


(5.18) 


(5.19) 


It  might  be  of  interest  to  note  here  that  limiting  forms  of  (5.19) 
lead  to  normalized  moments  of  either  the  lognormal  distribution  (r  -*•  °°) 
or  the  negative  exponential  distribution  (r  -*•  0),  i.e., 


£>  Jxn(n'1)/2  as  r  ~ 
£>n  (  n!  as  r  +  0, 


in  agreement  with  previous  theories.  (As  r  -*  0,  we  also  have  m  -*•  1  in 
(5.17)  to  get  this  result.) 

C.  Comnarison  of  Moments  with  Data: 


The  normalized  moments  (5.19)  depend  upon  three  parameters: 

(a)  .  the  reciprocal  of  the  normalized  variance  of  the 

m-distribution  denoted  by  m, 

(b)  .  the  normalized  second  moment  X  of  the  lognormal 

distribution,  and 

(c)  .  the  power  ratio  r  of  mean  intensities  of  the  specular 

to  diffuse  components. 

The  three  parameters  m,  X,  r  are  chosen  in  such  a  way  that  the  first  two 
normalized  moments  of  I  match  the  data  identically  and  the  third  moment 
is  matched  by  some  "best  fit"  criterion. 

Since  only  the  first  five  normalized  moments  were  actually  measured 
in  the  field,  we  separately  list  the  corresponding  theoretical  moments 
as  determined  by  (5.19), 


(5.20a 


P2  •+  4r  +  &2r 
(1  +  r)2 


(5.20b 


<I3>  _  U-  +  9y„r  +  9a„r2  +  a„r3 


(5.20d) 


and 

Pk+1  *  ^  +  k^vk»  k  "  1»2,3,4, 


this  last  result  being  derived  from  (5.17). 

The  measured  normalized  moments  <In>/<I>n  for  n  *  3,4,5  as  a  function 

2  2 

of  the  normalized  second  moment  <1  >/<I>  are  shown  in  Fig.  5-3  along 
with  theoretical  curves  predicted  by  the  present  G-distribution  (Eq. 
(5.20a)  -  (5.20e)).  In  Fig.  4-8  the  same  data  was  displayed  along  with 
theoretical  curves  predicted  by  the  lognormal  distribution  and  K- 
distribution. 

Although  the  moments  as  given  by  Eq.  (5.19)  depend  upon  three 
parameters,  we  found  that  a  reduction  to  one  effective  parameter  could 
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i-0.396x2  +,  5. 
0.2375X2  +  1. 


902x  -  4.506, 
4625x  -  1 .875, 


oT  _  <  2 
Lnl 

ar  -  >  2 
Lnl 


2  2 

where  x  *  <1  >/<I>  .  That  is,  we  found  the  normalized  second  moment 
associated  with  the  specular  component  varying  almost  linearly  with  the 
total  intensity  second  moment  whereas  a  relationship  involving  the 
normalized  second  moment  associated  with  the  diffuse  component  required 
a  quadratic  relation  which  further  depended  upon  the  value  of  the  log- 
intensity  variance.  The  parameter  m  is  then  calculated  from  (5.22)  and 
the  relation 


m  -  1/(^2  ~  !)• 


The  one  remaining  parameter  r  can  therefore  be  expressed  as  a  function 
of  X,  pj  and  x  utilizing  (5.20b)  from  which  we  deduce  the  relation 


2  -  x  +  /(VU  +  X  -  4)x.  +  4  -  Xu, 


Using  the  empirical  relations  (5.21)  and  (5.22)  along  with  (5.23) 

2 

and  (5.24),  the  variation  of  the  intensity  a  and  the  variation  of  the 

2 

power  ratio  r  as  functions  of  Oj^j  or  a  are  shown  in  Fig.s 
5-4  and  5-5,  respectively. 

We  should  hasten  to  point  out  that  the  empirical  relations  (5.21) 
and  (5.22)  were  not  arrived  at  by  applying  some  regression  relation  in 
the  hopes  of  finding  the  best  fit.  Since  the  actual  amount  of  data 
available  at  this  time  is  somewhat  sparse,  we  employed  visual  fits 
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5-5  *  7116  Power  ratio  of  mean  intensity  of  specular  to  diffuse 
components  of  scattered  laser  beam. 


based  upon  the  data  that  we  obtained  as  well  as  us' 4g  data  that  was 
earlier  published  by  Parry  and  Pusey  [9]*  It  is  likely  that  some 
modification  of  these  relations  will  result  when  more  data  is  available 
so  that  a  better  fit  of  the  theoretical  curves  can  be  obtained. 


VI.  BINARY  UNION  DECODING 


The  Binary  Union  Decoder  is  a  decoding  technique  used  to  recognize 
block  codes  which  are  transmitted  through  a  fading  channel.  The  MILES 
system  uses  such  a  decoder  in  order  to  defeat  the  fading  of  a  laser 
beam  signal  caused  by  scintillation  (or  other  scattering  mechanisms). 

The  decoder  consists  of  one  storage  register  that  acts  like  a 
memory  for  a  word.  The  next  received  word  is  then  OR’ed' bit  by  bit 
with  the  word  that  preceded  it  and  which  was  stored  in  the  memory 
register;  The  result  is  a  word  which  is  then  decoded.  This  technique 
is  designed  to  make  up  for  the  loss  of  one  half  of  a  word  so  long  as 
the  same  half  of  the  word  is  not  loiit  in  each  word  transmission. 

A  conceptual  block  diagram  of  the  decoder  is  shown  in  Fig.  6-1, 
illustrating  that  each  of  the  bits  in  the  stored  word  is  OR'ed  with 
its  corresponding  bit  in  the  current  incoming  word.  The  fact  that 
each  bit  is  used  twice,  once  when  it  enters  the  receiver  and  once  after 
it  has-  subsequently  been  stored  in  memory,  means  that  the  OR'ed  words 
shown  schematically  in  the  right-hand  column  are  not  independent  of 
one  another.  It  is  precisely  this  dependence  between  the  OR'ed  words 
that  was  ignored  in  an  earlier  analysis  of  the  decoder  and  thereby 
lead  to  an  overly-optimistic  prediction  of  the  MILES  System  performance 
121]. 

In  this  section  of  the  report  we  present  a  detailed  analysis  of 
the  Binary  Union  Decoder  taking  into  account  the  dependence  between 
OR'ed  words.  We  consider  MILES  block  code  words  with  weight  W  and  N+l 
repeated  word  transmissions.  Special  cases  to  be  considered  are  the 
following: 
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Case  I:  All  bits  in.  each  word  have  equal  probability  p  and  the  word 
detection  probability  computed  for  N  =  4  and  N  =  8.  (This  case 
simulates  the  situation,  wherein  the  statistics  of  the  propagation 
path  of  the  laser  beam  remains  constant  during  the  entire  message 
transmission.)  These  results  are  compared  with  the  current  MILES 
analysis  [21]  where  independence  between  OR'ed  words  was  assumed. 

Case  II:  All  bits  in  the  first  half  of  each  word  have  equal  probability 
p  and  all  bits  in  the  second  half  of  the  word  have  equal  probability 
q,  where  q  =  3p/4.  (This  case  models  the  fading  of,  half  a  word 
length-,  and  since  the  MILES  word  is  3.63  milliseconds  in  duration, 
it  corresponds  to  the  case  of  1.8  milliseconds  average  fade  time.) 

Case  III:  All  bits  in  the  first  half  of  each  word  have  equal  probability, 
and  this  probability  is  the  same  as  associated  with  all  bits  in 
the  last  half  of  the  preceding  word.  Thus  half  the  word  has 
probability  p  and  the  other  half  has  probability  q  =  3p/4.  (This 
case  simulates!  a  fading  of  3.63  milliseconds  but  only  half  of 
each  word  would  be  affected). 

Case  IV:,  All  bits  in  the  first  one-third  of  each  word  have  equal 
probability  p,  all  bits  in  the  next  one-third  of  each  word  have 
equal  probability  q,  and  all  bits  in  the  last  one-third  of  each 
word  have  equal  probability  r.  The  relationship  between  the 
probabilities  was  chosen  as  q  =  3p/4  and  r  =  *jp.  (This  case 
simulates  a  fade  length  of  about  1  millisecond.) 


Finally,  we  make  a  comparison  of  the  Binary  Union  Decoder  with 
the  Simple  Word  Recognition  Decoder.  The  Word  Recognition  Decoder  gives 
an  output  signal  only  after  it  has  recognized  the  received  bit  pattern 
as  a  valid  word.  Again  we  consider  MILES  block  code  words  with  weight 
W  and  N+l  repeated  word  transmissions.  The  special  cases  considered  here 
are  the  following: 

Case  V:  All  bits  in  each  word  have  equal  probability  p. 

Case  VI:  All  bits  in  the  first  half  of  the  word  have  equal  probability 
p  and  all  bits  in  the  second  half  have  equal  probability  q,  where 
q  =  3p/4. 

Case  VII:  One-third  of  the  bits  of  each  code  word  has  probability  p, 
one-third  has  probability  q  and  one-third  has  probability  r.  We 
again  assume  q  =  3p/4  and  r  =  ^p. 
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(6.2) 


where  1  £  i  S  N  and  1  £  j  £  W.  The  letter  N  denotes  the  number  of  unions 
and  thus  N+l  is  the  number  of  repeated  words  transmitted,  and  the  letter 
W  is  the  weight  of  the  code  word  (i.e.,  the  number  of  "ones”).  How  if 
only  one  good  word  is  required  for  a  successful  decode  (i.e.,  hit), 
then  mathematically  we  wish  to  determine  the  probability  associated  with 
N  unions  of  the  "A’s", 

N 

P(AjUA2u  ...  uAn>  P(UA.),  (6.3) 

i=l 


where  each  Aj  denotes  a  word  of  weight  W.  By  considering  only  the  "ones" 
and  not  the  "zeros" ,  we  are  considering  only  the  fading  loss  and  not  the 
more  general  case  of  noise  which  could  change  a  "zero"  to  a  "one"  as 
well  as  a  "one"  to  a  "zero". 

If  we  first  treat  the  special  case  of  only  three  repeated  words 
being  transmitted,  then  we  will  be  interested  in ‘computing 

PtAjU  A2)  =  P(Aj)  +  P(A2)  -  P(A1HA2).  (6.4) 

In  terms  of  the  OR'ed  bits  a..,  the  intersection  in  (6.4)  is 

A,«A2=  (aun  aJ2n  ...  ««lW)  n  (cf21*'a22n  *•*  fta2\P»  (6^) 
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or  by  regrouping, 


With  the  above  special  case  serving  as  a  model,  we  can  now  gen¬ 
eralize  the  results  to  the  case  when  N+l  repeated  words  are  transmitted. 
Here  we  find  the  probability  of  successful  decoding  one  word  out  of 
N+l  received  is 

P(  U  A  )  =  S  -  S  +  S  -  ...  +  (-1)N+1S  (6.13) 

3*1  i  1  1  5  N 

where  the  S  terms  are 

Sj  P(A.)  =(*)[l-(l-p)(l-p)]W 
N 

S2  =  Ep(Ai«  A  )  =ffi[l-(l-p)(l-p2)]W 
i  i,j-l  3  W 


50 


s 

mv* 

L 

E 

i 

1 

055  Q60  0.65  070  Q; 

BIT  PROB 

6-2.  Probability  of  detection  curves  for  the  case  when  all  bits 
have  equal  probability  and  N  =  4  or  N  =  8.  The  upper  curve  for 
each  N  corresponds  to  the  previous  MILES  analysis  given  by 
Eq.  (6.15)  and  the  lower  curve  for  each  N  corresponds  to  the 
present  analysis  based  on  Eq.  (6.14). 


S3I1P1 


N 

S„*  E PC*!"  »Ak)  =(“)fl-«-p)(l-p")]W- 

!)•••} k* 1 
tf...*k  m  terms 

Hence  the  probability  we  seek  is  given  by 

P(  "i  Aj)  =  (-l)m+1  ^[1-(1-P)(l-Pm)]  W-  (6.14) 

The  above  approach  to  the  analysis  is  an  accurate  accounting  of  the 
Biary  Union  Decoder.  It  might  be  pointed  out  that  if  the  unioned  words, 
A^,  were  assumed  independent,  as  is  the  case  in  the  current  MILES 
analysis  [21],  the  probability  of  a  successful  decode  would  be 

P(  U  A  )  =  1  -  (l  -  [1  -  (l-p)2]W}N+1  ,  (6.15) 

i*l  1 

where  N+l  is  again  the  number  of  repeated  words  and  W  is  the  weight  of 
the  word.  A  comparison  of  the  probabilities  generated  by  Eq.'s  (6.14) 
and  (6.15)  is  shown  in  Fig.  6-2.  Here  we  see  that  the  probabilities 
resulting  from  Eq.  (6.15)  are  much  higher  than  corresponding  proba- 
lities  given  by  Eq.  (6.14),  which  are  a  direct  result  of  ignoring  the 
dependence  of  the  unioned  words  A^. 

B.  Case  II  -  Half  Word  Fade  Time 

If  we  assume  a.  .  has  probability  p  for  j  =  l...,W/2  and  proba- 
bility  q  for  j  =  W/2  +  1,...,W,  then  Eq.  (6.7)  becomes 
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Fig.  6-3.  Probability  of  detection  curves  (lower  two)  for  the  case 
when  the  first  half  word  has  probability  p  and  the  second  half 
has  probability  q,  and  N  *  4  or  N  *  8  (see  Eq.  (6,19)).  The 
upper  two  curves  correspond  to  the  special  case  when  all  bits 
have  equal  probability  p  as  given  by  Eq.  (6.14). 


PROB  OF  FADE 


(6,16) 


W/2  W 

p(vv  =  if  H-u”  »2j> 

rl  j-f1 

Using  the  result  (6.10),  we  find 

P(A1"  A2)  =  [1  -  (1-pKl-p)2)]  W/2  [l  -  (l-q)(l-q2)]  W/2»  (6.17) 

and  therefore  it  follows  that 

p(VV  =  S  (“1)m+1  (m)fC1  "  (i-PHl-P®)]  [1  "  (l-q)(l~qm)]f  W/2-  (6.18) 

Generalizing  Eq.  (6.18),  we  have 

N 

P(  U  At)  =  J^(-l)m+1  ^j[l-(l-p)(l-pm)][l-(l-q)(l-qm)]lW/2*  (6.19) 

Eq.  (6.19)  is  plotted  in  Fig.  6-3  for  the  special  case  when  the 
bit  probabilities  satisfy  q  =  3p/4  and  N  =  4  or  N  =  8.  Corresponding 
curves  for  the  case  when  all  bits  have  probability  p  (see  Eq.  (6.14)) 
are  illustrated  here  for  comparison  purposes.  The  relative  bit  proba¬ 
bilities  chosen  were  picked  only  as  typical  of  fading  conditions  and 
also  to  demonstrate  the  effect  of  the  Binary  Union  on  the  word  detection 
probability. 

C .  Case  III  -  Full  Word  Fade  Time  -  Half  Word  Lost : 

Here  we  assume  P(a.  .)  =  p,  j  =  l,...,W/2  and  P(a1  .)  =  q, 

j  =  W/2  +1,...,W.  Then  P(a9  •)  =  q,  j  =  l,...,W/2  and  P(a9  .)  =  p, 

'••J 

j  =  W/2  +1,...,W,  and  so  on  where  every  other  word  has  the  same  proba¬ 
bilities  assigned  to  each  bit. 
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and  likewise 

P(A2)  =  [1  -  (l-q)(l-p)]  W 

Also,  we  find  that 

W/2  W 

F(Ai"V  *  "  .  VP(oij  "  V 

3*1  3*2+1 

=  [2(p+q-pq)-(2p+q-2pq-p2+p2q)]  W^2  [>2(p+q-pq)-(2q+p-2pq_q2+pq^)]  W^2 
=  [1  -  (l-q)(l-p2)]W/2  [1  -  (l-p)(l-q2)]  W/2 


(6.23) 
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Fig.  6-4.  Probability  of  detection  curves  for  the  case  when  the  first 
half  of  every  other  word  has  probability  p  and  the  second  half 
has  probability  q,  alternated  with  those  words  where  the  first 
half  has  probability  q  and  the  second  half  has  probability  p,  and 
N  «  4  or  N  =  8  (see  Eq.  (6.25)).  The  dashed  curve  for  each  N 
corresponds  to  die  probability  of  detection  when  each  word  has 
the  same  p,q  probability  assignment  as  the  previous  word  (see 
Eq.  (6.19)). 
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6-5.  Probability  of  detection  curves  for  the  case  when  the  first 
one- third  of  the  word  has  probability  p,  the  next  one-third  has 
probability  q,  and  the  last  one-third  has  probability  r,  where 
q*3p/4  and  r*^p,  and  N  =  4  or  N  =  8  (see  Eq.  (6.26)).  The  upper 
two  curves  correspond  to  the  special  case  when  all  bits  have 
equal  probability  p  as  given  by  Eq.  (6.14)>. 
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Therefore,  Eq.  (6.20)  now  takes  the  form 


i 

f 


u 


u 

U 


I 


MW 


P(A1  V  ^+1  (‘jilH1  W!  [i-(i-p)U-qm3>  w/2,  (6.24) 

and  generalizing  to  N+l  repeated  words  yields 
N  N  /  v 

p(iUA.)  =^(-l)mKi^jl[l“(l-q)(l-pm3  [l-(l-p)  (l-qm3  J  W/2.  (6.25) 

Probability  of  detection  curves  using  Eq.  (6.25)  for  the  special 
case  when  q  =  3p/4  and  N  =  4  or  N=  8'  are  shown  in  Fig.  6-4.  The  upper 
curves  for  each  N  illustrate  similar  probabilities  of  detection  when 
each  word  has  the  same  probabilities  p  and  q  as  determined  by  Eq.  (6.19), 
rather  than  every  other  word  as  is  the  case  of  present  interest* 


D.  Case  IV  -  Third  of  a  Word  Faded; 

Finally  we  consider  when  P(a..  .)  =  p,  j  =  l,...,W/3,  P(a.  .)  =  q, 

i»J  X»J 

j  =  W/3+1, . . . ,2W/3,  and  P(a,  .)  =  r  for  j  =  2W/3+1 , . . . ,W.  Following  the 
same  procedure  as  outlined  in  Section  B,  we  find  that 

p ( U Ai )  1 )m+ 1 (m jV1- ( 1 "P) ( 1"Pm^  ll-(l-q)(l-qm1  [l-(l-r)(l-rra)]|W/3‘  (6.26) 

m=l  '  ' 

Probability  of  detection  curves  for  this  case  are  illustrated  in 
Fig.  6-5  for  the  special  case  when  q  =  3p/4  and  r  =  *jp.  Again  these 
values  represent  typical  probabilities. 

E.  Case  V  -  All  Bits  Equal  Probability: 

If  we  let  the  probability  of  detection  of  a  single  bit  be  p,  the 
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word  detection  probability  for  a  word  of  weight  W  is 


p  w 

word  ~  ^ 

The  probability  of  missing  a  word  is 


(6.27) 


(6.2S) 


and  the  probability  of  missing  every  one  of  N+l  repeatedly  transmitted 
words  is 


(1  -  PW)N+1. 


(6.30) 


F.  Case  VI  -  One  Half  Word  Fade: 

Suppose  p  is  the  probability  of  detecting  half  of  the  bits  in  a 
word  and  q  is  the  probability  of  detecting  the  second  half  of  the  bits 
in  the  same  word.  The  probability  of  detecting  the  word  is 


P 


word 


(6.31) 


and  hence  the  probability  of  detecting  one  word  in  N+l  repeated  trans¬ 


missions  is 


G.  Case  VII  -  One-Third  of  a  Word  Faded: 

Let  the  probability  of  each  third  of  a  word  be  p,  q,  and  r, 
respectively.  Following  the  analysis  in  Section  F,  we  deduce  that  the 
probability  of  detecting  a  single  word  when  N+l  are  transmitted  is 

PD  =  1  -  [1  -  (pqr)W/3]  N+1.  (6.33) 

Fig.  6-6  shows  the  plots  of  Eq.'s  (6.30),  (6.32)  and'  (6.33)  where 
we  have  chosen  q  =  3p/4  and  r  =!*p.  These  choices  represent  typical 
fading  depths  and  rates  as  discussed  earlier.  Comparison  of  the  results 
here  for  the  Word  Recognition  Decoder  and  those  in  previous  sections 
for  the  Binary  Union  Decoder  clearly  reveals  the  effectiveness  of  the 
Binary  Union  Decoder.  Observe,  however,  that  the  Binary  Union  Decoder 
is  most  effective  by  comparison  with  the  Word  Recognition  Decoder  when 
all  bits  have  equal  probability. 
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NON*  BINARY  UNION 


•  -  Alt  bits  equal  prob.  p 

©  •  Bit  prob.  =  p,q  where  q=  -f-p 

0  -  Bit  prob.  =p,q,r  where  q=  ,  -2-p  and  r  = 


Z  3- 

G 

6  *' 
£  *■ 
a  «■ 

u.  8' 

O  *. 


GO  9 

9 

‘  i 
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Fig.  6-6.  Probability  of  detection' curves  for  the  Word  Recognition 
Decoder  corresponding  to  either  N  *  4  or  N  =  8.  The  upper  two 
curves  correspond  to  Eq.  (6.30)  where  all  bits  have  equal  proba¬ 
bility,  the  middle  two  curves  correspond  to  that  case  when  half 
the  received  word  has  probability  p  and  half  has  probability  q 
(Eq.  (6.32)),  and  the  lower  two  curves  correspond  to  that  case 
when  the  word  probability  is  split  into  thirds  (Eq.  (6.33)). 
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VII.  DISCUSSION 

The  purpose  of  this  study  was  to  evaluate  the  potential  for  using 
MILES  equipment  for  long  range  marksmanship  gunnery  training.  The 
current  MILES  equipment  was  designed  as  a  tactical  engagement  laser 
weapon  fire  simulator.  In  the  development  of  the  system  detailed 
analysis  was  performed  to  predict  the  optical  communication  process. 
The  Biliary  Union  Decoder  was  designed  to  increase  the  detection  proba¬ 
bility  of  received  MILES  words  under  very  specific  fading  conditions. 
That  analysis  was  performed  assuming  short  range  propagation  path  and 
weak  scattering  in  the  path.  Also,  simplifications  were  made  in  the 
analysis  of  the  MILES  decoder  circuit  which  are  valid'  for  the  short 
range  problem  but  which  proved  unreliable  and  misleading  for  the  long 
range  problem.  That  is,  if  the  short  range  (200  meters)  analysis  is 
extended  to  the  longer  ranges  the  probability  of  detection  curves  that 
result  are  overly  optimistic. 

In  order  to  definitively  evaluate  the  MILES  potential  for  long 
range  marksmanship  gunnery  training,  the  following  tasks  are  required: 

1.  Evaluate  currently  available  mathematical  models  and  alter 
or  create  new  models  for  long  range  optical  propagation. 

2.  Experimentally  evaluate  MILES  long  range  detection  proba¬ 
bility  (kill  probability). 

3.  Develop  appropriate  and  necessary  add-ons  for  MILES 
equipment  which  will  allow  it  to  accurately  simulate  long 
range  gun  firing  for  marksmanship  training. 

A.  Completed  Modelling  Tasks: 

During  the  current  study  we  have  completed  the  following  aspects 


of  the  mathematical  modelling: 


(a) .  Experimentally  determined,  in  well  documented  field 

experiments,  conducted  at  NASA's  Kennedy  Space  Center, 
that  the  Lognormal  Statistical  Model  for  scintillation 
scattering  is  valid  only  for  ranges  less  than  200  meters 
(or  weak  turbulence  conditions). 

(b) .  The  field  experiments  also  proved  that  the  variance  cf  the 

statistically  fluctuating  optical  signal  at  long  ranges 
does  not  saturate  at  the  low  values  (0.7)  used  in  the  MILES 
analysis.  The  actual  variance  at  long  ranges  may  be  as 
high  as  4  dr  5. 

(c) .  The  K-distribution  was  validated  as  a  reasonable ^vnodel  for 

long  ranges  although  it  is  not  theoretically  applicable 
for  short  ranges  when  the  normalized  variance  is  less  than 
unity. 

(d) .  The  cumulative  fading  (or  cumulative  distribution)  was 

computed  for  the  K-distribution  and  for  the  lognormal 
distribution  using  the  same  variance  in  both  cases  (see 
Fig.  7-1).  The  lognormal  model  predicts  lower  fading 
at  long  ranges  as  compared  with  the  K-distribution  leading 
to  a  more  optimistic  view  of  system  performance. 

(e) .  A  universal  statistical  model,  called  the  G-distribution, 

was  developed  that  would  accurately  predict  scattering 
effects  at  both  short  range  and  long  range.  Due  to 
complexities  of  the  mathematical  functional  form  of  the 
G-distribution,  only  the  statistical  moments  have  been 
computed  with  the  distribution  at  this  time.  Further 
analysis  and  refinement  of  this  distribution  is  required 
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to  cast  it  into  a  form  that  will  be  propitious  to  the 
computation  of  cumulative  fading. 

In  addition  to  the  mathematical  modelling  of  the  intensity  fluc¬ 
tuations  of  the  optical  beam,  an  analysis  was  performed  on  the  Binary 
Union  Decoder  employed  in  the  MILES  system.  Although  the  current 
analysis  was  reasonably  accuracy  for  short  ranges  [21],  it  leads  to 
a  considerably  over-optimistic  performance  of  the  MILES  system  for 
longer  ranges.  The  analysis  performed  in  this  study  for  the  decoder 
made  use  of  the  lognormal  distribution  as  well  as  the  K-distribution 

to  illustrate  the  differences  in  the  predicted  performance  of  MILES 

/ 

at  long  ranges  for  both  statistical  models  (see  Fig.  7-2). 


B.  Evaluation  of  MILES  Equipment: 

As  of  the  completion  of  this  portion  of  the  study,  there  has  been 


no  MILES  equipment  available  for  testing.  If  and  when  the  equipment 
becomes  available,  it  will  be  taken  to  NASA's  Space  Shuttle  Runway  for 


the  purpose  of  conducting  experiments  in  the  hopes  of  determining 


detection  (kill)  probability  zones  as  a  function  of  range. 


MILES  Equipment  Add-On: 


The  major  portion  of  the  marksmanship  study  concerning  MILES 
hardware  add-ons  was  conducted  by  Code  N-73  of  NTEC's  research  divi¬ 
sion.  The  UCF  research  team  worked  with  Code  N-73  in  the  early  stages 
of  the  study  to  help  develop  add-on  hardware  concepts  that  would  be 
MILES  compatible  and  allow  MILES  to  be  used  in  marksmanship  gunnery 
training.  The  concepts  that  evolved  from  the  preliminary  discussions 
were: 
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1.  Phase  locked  loop  (PLL)  receiver. 

2.  Digital  filter  (DF)  receiver. 

3.  Continuous  wave  (CW)  laser,  transmitter  pulsed  at  high  rate. 

Since  the  MILES  transmitters  operate  at  a  basic  transmission  rate 

of  3  khz,  the  high  pulse  rate  CW  laser  could  be  configured  so  as  to 
transmit  20  pulses  during  the  slot  time  for  one  MILES  bit  pulse.  The 
receivers  would  then  lock  onto  the  20  pulses  and  output  one  pulse  at 
the  MILES  bit  pulse  period.  Hence  the  proposed  transmitter  and 
receiver  would  appear  as  standard  MILES  transmitter  and  receiver  to 
the  present  MILES  encoder  and  decoder  electronic  systems.  In  order  to 
drive  the  transmitter  and  receiver  developed  by  Code  N-73  of  NTEC, 
the  UCF  research  team  developed  the  software  so  that  the  microprocessor 
system  PROMPT  80  could  be  used  as  MILES  encoder.  The  MILES  decoder 
system  can  be  simulated  using  the  MINC-11  computer.  Appendix  C  lists 
the  microprocessor  program  for  simulating  the  MILES  encoder. 

The  above  techniques  were  being  developed  simultaneously  and  in 
conjunction  with  the  UCF  study  of  the  mathematical  modelling.  These 
techniques  are  well  known  in  communication  theory  and  appear  to  be 
very  promising  as  techniques  to  allow  MILES  to  be  used  for  marksman¬ 
ship  gunnery  training. 


APPENDIX  A 


The  Characteristic  Function 

Let  us  write  the  time-independent  portion  of  the  field  as 

E  -  Ae19  +  Re1^  -  Ex  +  E2-  (A-l) 

Since  E^  and  are  complex  functions,  we  define  their  characteristic 
functions  as  two-dimensional  Fourier  transforms.  For  E^,  this  becomes 

2  it 

C  (u,v)  -/  /eluAcos9  +  iv^8in9p(A,0)  dGdA,  (A-2) 

o  o 

where  p(A *0)  is  the  joint  density  function  of  A  and  6.  If  we  assume 
that  A  and  0  are  independent,  and  further  that  6  is  uniformly  dis- 
tirbuted,  then 

p(A,0)  *  p(A)p(0)  *  2^"  P<A) 

and  (A-2)  can  be  written  in  the  form 

00  2tt 

c^(u,v)  »  /  p(A)  j  eiAzcos^9  "  do  dA,  (A-3) 

o  o 

where 


and 

a  =  arctan(v/u). 

Using  the  integral  formula  ([17.1,  p.  360) 

2tt 

i-  jreixcc»rf  u  .  J'(x)  (A-4) 

0 

it  becomes  clear  that  depends  only  upon  the  radial  variable  z,  i.e., 

00 

C1<Z)  -  j  p(A)J0(Az)  dA,  (A-5) 
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r 


which  we  recognize  as  being  equivalent  to  the  Hankel  transform  of  p(A)/A. 
Under  similar  assumptions  for  R  and  the  characteristic  function  for 


leads  to 


C2(2) 


/  p(R)Jo<Rz)  dR  , 


(A-6) 


and  therefore  the  characteristic  function  for  E  -  can  be  expressed 


in  the  product  form 

00  00 

C(z)  -  C1(z)C2(z)  -  J  p(A)Jo(Az)dA  f  p(R)JQ(Rz)dR.  (A-7) 


The  pdf  for  E  is  defined  as  the  two-dimensional  Fourier  transform 
of  C(z)  given  by 


P(E)  '£// 


-iuX  -  ivY  a,  s  ,  , 

e  C(z)  du  dv> 


(A-8) 


where 


X  ■  Re(E)  *  A  cos  0  +  R  cos  6 


and 


Y  *  Im(E)-  ■  A  sin  0  +  R  sin 

Since  the  characteristic  function  is  formulated  in  terms  of  the  radial  variable 
z  it  is  convenient  to  write  the  exponential  function  in  (A-8)  in  terms  of 
polar  coordinates.  Thus  we  write 


uX  +  vY  «  |E|z  cos(^'  -  o) 


where 


|E| 


2  2 
X  +  Y 


and 


iJj  =  arctan(Y/X). 


Converting  the  double  integral  in  (A-8)  to  polar  coordinates  (z,a)  it  follows 
that 


p(E) 


.  1 


4tt 


f  f  zCMe'llElz  c°sW  "  a>  d«d« 


(A-9) 


0  0 


which  in  view  of  (A-4)  further  reduces  to 
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APPENDIX  B 

The  Meijer  G-Function 


In  Section  VA  we  developed  the  expression 

Pd)  -I  r(k+1)  cf3(fi|  1 


k-0  (kl)  rOn) 


1-1,  m,  k+l/ 


(B-l) 


as  the  pdf  associated  with  the  intensity  fluctuations  of  an  optical  beam 

21 

propagating  through  clear-air  turbulence.  By  »  we  mean  one  member  of 

the  generalized  function  known  as  the  Meijer  G-f unction  [19].  It  is  our 

intent  here  to  relate  this  function  to  more  standard  functions  used  in 

21 

engineering  applications.  In  particular,  we  will  show  that  G23^ 
related  to  the  generalized  Laguerre  functions. 

Using  the  relations  ([19],  vol.  1,  p.  209  and  p.  216) 


G23  x 


b,  c,  a 


•  -g  * 


(B-2) 


_20  /  \  ’-s(b+c-l)  -Hpc  ..  ,  . 

“I  b,  cj  6  Wb-c)<:,)’ 


(B-3) 


we  can  relate  a  particular  G23^*)  function  to  the  Whittaker  functions 

W  (.).  Then  employing  the  relation  ([19].  vol.  2,  p.  432) 
n ,  v 

WU-fn+Js,  ±u(x)  "  C-l)n  n!  L<2u)  (x)  e  */z,  (B-4) 


where  (.)  is  the  generalized  Laguerre  function,  we  have  one  special 
case  of  interest 


U  3  )  -  (-l)b'a  (b  -  .)!  xc  e-x  L<c-b>  <x). 

V  b,  c,  a  /  b'a 


(B— 5) 
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The  G-function  in  (B-l)  is  more  general  than  that  in  (B-5) .  There¬ 
fore  in  order  to  utilize  (B-5)  we  need  to  develop  an  appropriate  recurrence- 
type  formula  which  relates  the  more  general  G-function  to  those  of  the  specific 
fora  in  (B-5).  Such  a  recurrence  formula  can  be  derived  by  making  use  of  the 
relation 

G23  (X  1  ^  1  )  ~  (k+j-n)G^/x|  ^  1  (xj  )  (B~6) 

\  1  k+j,  m,  k+1  /  \  I  k+j,  m,  k+1  /  \  »k+j+l,m,k+l/ 

•Which  is  a: .special  case  of  Eq.  (11),  p.  209,  vol.  1,  [19).  Starting  with 
jv ■  1  and  n  *  1,  repeated  application  of  (B-6)  yields 


k(k-l)G23  (x 


k+1, m, k+1 


)'2kG^(X  k+2,»,k+i) 

(.  3>1  )  • 

A  k+3,m,k+l/ 


while  in  general  it  can  be  shown  that 


„2I\  / 

°23  (x 


k+1,  m,  k+1  /  j 


i  <-!>*%*)*  i.  k+1'1 

5=0  'J/  •  \  1 2k-j+l,m,l 


.  (B-7) 


Recalling  (B-5) ,  this  last  result  leads  to 


r21/ v 
23\ 


\«  (-l)kxVxx(j)  j!(2k-j)lL^k+j"1)(x),  (B-8) 

a,  k+1/  j-0  4K‘3 


and  the  substitution  of  this  expression  with  x  *  ml/b  into  (B-l)  gives  us 

P(I)  -  e-I/b  I  I  (aei>L<A!i  j-k  L(»-2^-i)(Bl/b).  (B.9) 

r(m)bm  k=0  j-0  j![(k-j)!'J  3 


and  by  setting  lc  -  j  *  n, 


,  k 


j-0  jl  C(k-j) 


k-J)!r  d  kl  n“0  '  '  1 


n 


(-ir 

bk  k! 


bk(I/b), 


where  1^(0  is  the  kth  Laguerre  polynomial.  Thus,  by  setting  m  *  1 
in  (B-9)  and  using  the  above  relation,  we  see  that  (B-9)  reduces  to 


P(I) 


k«0 


(~l)k<A2k> 
bk  k! 


Lk(I/b) 


(B-10) 


2  2 

which  is  in  agreement  with  Eq.  (19)  in  C83  for  I  ■  E  and  b  »  2a^, 

That  is  to  say,  (B-10)  is  the  pdf  that  would  result  under  the  assumption 
that  the  diffuse  component (s)  have  an  amplitude  R  that  is  Rayleigh 
distributed  rather  than  the  more  general  m-distribution  as  assumed  for 
our  model  leading  to  (B-9) . 


APPENDIX  C 


Mircoprocessor  Simulator 
and  High  Data  Rate  Receiver 
C.l  Transmitter  and  Encoder: 

In  order  to  conduct  experiments  on  the  high  data  concepts  for 
MILES  add-on's  it  is  required  that  the  MILES  encoder-transmitter  be 
simulated.  The  encoder  transmitter  design  required  that  the  system 
be  as  flexible  as  possible.  This  design  had  to  allow  for  a  variety 
of  experiments  on  different  types  of  modulation  and  block  coding 
schemes.  Because  of  the  diversity  of  digital  block  codes,  it  was 
decided  that  a  general  purpose  microprocessor  computer  should  be  used 
to  generate  these  codes.  The  different  proposed  concepts  also  required 
that  the  transmitter  be  able  to  transmit  block  codes  either  by  Frequency 
Shift  Keying  (FSK)  or  by  On-Off  Keying  (00K)  modulation  formats. 

The  use  of  a  computer  for  generating  block  codes  requires  the 
computer  to  interface  with  a  teletype  and  to  have  an  I/O  port  that  is 
easily  accessible.  The  computer  used  for  this  task  was  the  Prompt 
80,  an  Intel  8080  (CPU)  computer,  which  is  capable  of  driving  the 
necessary  teletype  and  has  an  easily  accessible  I/O  port.  The  Prompt  80 
is  limited  to  IK  of  R.A.M.  memory,  but  because  it  doesn't  have  an  on-board 
assembler,  programs  must  be  hand-assembled  and  keyed  in  by  hand  via  the 
front  panel  keyboard. 

The  code  generation  program  had  to  produce  block  codes  that  were 
MILES  compatible.  The  MILES  block  codes  consist  of  32  different  11-bit 
words.  These  words  are  sent  at  a  bit  rate  of  333  sec  per  bit.  Each 
word  is  then  repeated  and  re-transmitted  a  specified  number  of  times. 

For  the  Dragon  missle,  the  Tow  missle  and  the  Viper,  a  special  subroutine 


was  needed.  To  reduce  the  number  of  program  subroutines  required, 

the  number  of  repeats  of  the  word  transmission  was  entered  via  the 

teletype.  A  Print  subroutine,  an  ASCII  binary- to-binary  subroutine, 

an  ASCII  decimal-to-binary  subroutine,  and  an  .output  subroutine  were 

the  major  subroutines  written.  The  Print  subroutine  printed  the  messages 

on  the  teletype.  The  ASCII  binary-to-binary  subroutine  and  the  ASCII 

decimal-to-binary  subroutine  were  two  required  conversion  routines.  The 

ASCII  binary-to-binary  subroutine  was  used  to  enter  the  11-bit  code  word 

while  the  ASCII  decimal-to-binary  subroutine  was  used  to  enter  the  number 

* 

of  times  the  code  would  be  repeated.  The  maximum  repetion  of  the  code 
was  set  at  99.  The  output  subroutine  serially  clocked  the  11-bit  code 
words  bit  by  bit  at  a  rate  of  333  p  sec  per  bit.  The  program  listing 
is  given  in  Table  C-l. 

The  output  of  the  computer  was  used  to  control  the  modulator 
which  in  turn  was  used  for  the  two  different  modulating  formats  (FSK 
and  00K).  The  circuit  developed  allowed  for  both  methods  to  be  utilized 
with  a  minimum  number  of  components  (see  Fig.  C-l),  and  switches  from 
FSK  to  OOK  by  the  throw  of  a  switch.  The  modulator  circuit  performs 
in  the  following  way.  The  flip-flop  and  nand  gates  1  and  2  are  used 
to  divide  the  input  carrier’s  frequency  by  two.  These  three  logic 
devices  allow  the  carrier's  frequency  to  be  divided  by  two  and  yet  be 
independent  of  the  carrier's  pulse  width  (minimum  of  60  nsec  pulse 
width).  With  the  switch  closed,  there  is  no  output  from  nand  gate 
2  and  hence  its  output  remains  low  (OV).  The  computer  code  controls 
which  carrier  signal  is  allowed  at  the  output  of  nand  gate  6;  this  is 
accomplished  with  nand  gates  J  and  5.  If  the  computer  output  is  low 
(OV) ,  a  high  (+5V)  appears  at  the  input  to  nand  gate  3  and  a  low  (OV) 
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Fig.  C-l(a).  Modulator  configuration  for  FSK  and  OOK  modulation  of  MILES 
codes. 
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Fig.  C-l(b),  Timing  diagram  of  modulator. 
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appears  at  the  input  to  nand  gate  5.  A  high  (+5V)  at  the  input  to 
aand  gate  3  allows  the  carrier  (with  frequency,  divided  by  two)  to 
be  present  at  the  output.  If  the  computer  code  is  high  (+5V), 
the  input  to  nand  gate  3  is  low  (OV)  and  the  input  to  nand  gate  5 
is  high  (+5V).  Again  this  peraits-  the  carrier  (with  frequency  divided 
by  one)  to  be  present  at  the  output.  The  result  is  FSK  modulation. 

If  the  switch  is  closed,  the  circuit  operates  as  described  above 
except  when  the  computer  code  is  low  (OV) .  With  the  switch, closed 
there  is  no  output  from  the  modulator  when  the  computer  output  is 
low  (OV)  which  Results  in  00K  modulation.  The  output  of  the  modulator 
.is  tihen  fed  to  the  driver  circuitry  which  drives  the  laser  diode.  The 
laser  is  then  pulsed  at  the  carrier  frequency  or  carrier  frequency 
divided  by  two  depending  upon  the  state  of  the  computer  output  and  the 
modulator  switch  position. 


TABLE  C-l 


MICROPROCESSOR  8080  PROGRAM  FOR  CONTROLLING 
TRANSMITTER  TO  TRANSMIT  MILES  CODE 


ADDRESS 

ASSEMBLY 

OP  CODE 

3D00 

Start  , 

LXI  SP,  //3CF0 

31F03C 

3D03 

LXI  HL,  //ME SI 

21D03E 

3D06 

CALL  OUT 

C0003E 

3D09 

MVI  d,  m 

1609 

3D08 

MVL'B,  #0 

0600 

3D0D 

Start  1: 

DCR  D 

15 

3D0E. 

JZ  Start  2 

CA773D. 

3D11 

CALL  1BC2 

CD2B3E 

3D14 

JMP  Start  1 

C30D3D 

3D17 

Start  2: 

MOV  E,B 

58 

3D18 

MVI  B,  //O 

0600 

3D1A 

MVI  D,  #4 

1604 

3D1C 

Start  3: 

DCR  D 

15 

3D  ID 

JZ  Start  4. 

CA263D 

3D20 

CALL  ilBC2 

CD2B3E 

3D23 

JMP  Start  3 

C31C3D 

3D26 

Start  4: 

LXI  HL,  //Temp  1 

21F13C 

3D29 

MOV  M,E 

73 

3D2A 

I NX  H,L 

23 

3D2B 

MOV  M,B 

70 

3D2C 

Start  5: 

LXI  HL,  //MES2 

21033F  . 

3D2F 

CALL  OUT 

CD003E 

78 


ADDRESS 

ASSEMBLY 

OP  CODE 

3D32 

MVl  E,  00 

1E00 

3D34 

CALL  DECB 

CD333E 

3D37 

LXI  #Temp  3 

21F33C 

3D3A 

MOV  M,E 

73 

3D3B 

Start  6: 

LXI  HL,  MES3 

21313F 

3D3E 

CALL  OUT 

CD003E 

3D41 

CALL  IN 

CD0C3E 

3D44 

CPI  #53 

FE53 

3D46 

JNZ  Start  6 

C2003D 

3D49 

CALL  SHIFT 

CD7D3E 

3D4C 

CALL  OUT  1 

CD593E 

3D4F 

LXI  HL,  #MES4 

21543F 

3D52 

CALL- OUT 

CD003E 

3D55 

JMP  Start 

C3003I 

3E00 

OUT: 

MOV  A,M 

7E 

3E01 

CPI  #00 

FEOO 

3E03 

RZ 

C8 

3E04 

MOV  C,A 

4F 

3E05 

CALL  07FA 

CDFAO 

3E08 

I  NX  HL 

23 

3E09 

JMP  OUT 

C3003 

3E0C 

IN: 

CALL  0729 

CD290 

3E0F 

MOV  C,A 

4F 

3E10 

CALL  07FA 

CDFAO 

3E13 

ANI  #7F 

E67F 

79 


ADDRESS 

ASSEMBLY 

OP  CODE 

3E15 

RET 

C9 

3E16 

IBC: 

CALL  IN 

CD0C3E 

3E19 

cpi  m 

FE49 

3E1B 

JZ  #MISS 

CA003D 

3E1E 

CPI  #32 

FE32 

3E20 

JNC  #Start 

D2003D 

3E23 

CPI  #30 

FE30 

3E25 

JC  #Start 

DA003D 

3E28 

ANI  #1 

»  E601 

3E2A 

RET 

C9 

3E2B 

IBC2: 

CALL  IBC 

CD163E 

3E2E 

RAR 

IF 

3E2F 

MOV  A,B 

78 

3E30 

RAL 

17' 

3E31 

MOV  B,A 

47 

3E32 

RET 

C9 

3E33 

DECB: 

MV I  D,  #2 

1602  • 

3E35 

PUSH  DE 

D5 

3E36 

LXI,  D',0 

110000' 

3E39 

LOOP: 

CALL  IN 

CD0C3E 

3E3C 

CPI  #3A 

FE3A 

3E3E 

JNC  Start  5 

D22C3D 

3E41 

CPI  #30 

FE30 

3E43 

JC  Start  5 

DA2C3D 

3E46 

ANI  #0F 

E60F 

80 


ADDRESS 

ASSEMBLY 

OP  CODE 

3E48 

ADDE 

83 

3E49 

NOV  E,A 

5F 

3E4A 

POP  B,C 

Cl 

3E4B 

DCR  B 

05 

3E4C 

RZ 

C8 

3E4D 

PUSH  B 

C5 

3E4E 

XCHG 

EB 

3E4F 

DAD  H,L 

29 

3E50 

PUSH  H,L 

E5 

3E51 

DAD  H,L 

29 

3E52 

DAD  H,L 

29 

3E53 

POP  B,C 

Cl 

3E54 

DAD  B,C 

09 

3E55. 

XCHG 

EB 

3E56 

JMP  LOOP 

C3393E 

3E59 

OUT  1: 

LX1  HL,  Temp  1 

21F33C 

3E5C 

’mov  c,m 

4E 

3E5D 

LHLD  Temp  1 

2AF13C 

3E60 

XCHG 

EB 

3E61 

OUT  2: 

MOV  H,E 

63 

3E62 

MOV  L,D 

6A 

3E63 

MVI  B,  //B 

060B 

3E65 

MOV  A,H 

7C 

3E66 

ANI  #80 

E680 

3E68 

RAL 

17 

81 


ADDRESS 


ASSEMBLY 


3E69 

3E6A 

3E6C 

3E6D 

3E70 

3E72 

3E73 

3E76 

3E78 

3E79 

3E7C 

3E7D 

3E80 

3E81 

3E82 

3E83 

3E84 

3E85 

3E86 

3E87 

3E89 

3E8A 

3E8D 

3E83 

3E8F 


SHIFT: 


DELAY: 
DELAY  1: 


RAL 
OUT  E8 
DAD  H,L 
CALL  DELAY 
MVI  A, 00- 
DCR  B 
JNZ  OUT  3 
OUT  E8 
DCR  C 
JNZ  OUT  2 
RET 

LX1  HL,  Temp  2 

MOV  A,M 

RAR 

RAR 

RAR 

RAF" 

MOV  M,A 
RET 

mvi  a,  m 

DCR  A 
JNZ  DELAY 
INC  A 
INC  A 

MVI  A,  #14 


OP  CODE 
17 

D3E8 

29 

CD873E 

3E00 


s 

? 


2 

Jf 


ADDRESS 


ASSEMBLY 


OP  CODE 


3ED0 

Message  1: 

/CR  LF  Enter  11  Bit  Character 

Code  or  M  For  Missiles?  CR  LF/ 

3F03 

Message  2: 

/CR  LF  Enter  The  Number  of  Times 
To,  Repeat  Code?  CR  LF/ 

3F31 

Message  3: 

/CR  LF  Press  S  To  Start  Trans¬ 
mission.  CR  LF/ 

3FS4 

Message  4: 

/CR  LF  Code  Has  Been  Transmitted. 
CR  LF/ 

TABLE  C-2 
SYMBOL  TABLE 


MES  1 

ZZ 

3ED0 

TEMP  1 

ZZ 

3CF1 

MES  2 

zz 

3F03 

TEMP  2 

= 

3CF2 

MES  3 

— 

3F31 

TEMP  3 

s 

3CF3 

MES  4 

3F54 

SHIFT 

= 

3E7D 

START 

3D00 

IBC 

= 

3E16 

START  1 

3D0D 

IBC  2 

= 

3E2B 

START  2 

= 

3D17 

DECB 

3E33 

START  3 

3D1C 

LOOP 

3E39 

START  4 

ZZ 

3D26 

OUT  1 

3E59 

START  5 

3D2C 

OUT  2 

= 

3E61 

START  6 

= 

3D3B 

OUT  3 

3E65 

OUT 

=: 

3E00 

DELAY 

= 

3E87 

DELAY  1 

= 

3E88 

NOTE  1:  0729  is  the -address,  for  the  Prompt  80' s  input  subroutine, 

(input  one  ASCII  character) 

NOTE  2:  07FA  is  the  address  for  the  Prompt  80' s  output  subroutine, 
(output  one  ASCII  character) 
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Fig.  C-2.  Circuit  connections  for  PLL  chip. 


PLL  INPUT 


Fig.  C-3.  Output  waveforms  when  PLL  is  locked  onto  signal  and  when  it 
is  free-running. 
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In  order  to  check  the  computer  program  the  Tektronix  7704  Digital 
Processing  Oscilloscope  and  the  Hewlett  Packard  5333A  frequency  counter 
were  used.  To  verify  the  program  a  code  word  of  10000000000  was  entered 
and  repeated  indefinitely,  which  corresponded  to  a  frequency  of  273  Hz. 

The  first  bit  was  high  for  332  p  sec  leading  to  a  corresponding  word 
frequency  of  271.1  Hz.  This  meant  there  was  an  error  of  0.6%  between 
the  MILES  and  Prompt  80  code  generators. 

C.2  Receivers: 

A.  Phase  Locked  Loop  (PPL): 

A  phase  locked  loop  is  an  electronic  feedback  serve  loop  con¬ 
sisting  of  a  phase  detector,  a  low  pass  filter  and  a  voltage  controlled 
oscillator.  The  controlled  phase  of  the  oscillator  makes  this  system 
capable  of  tracking  and  locking  onto  a  .received  periodic  signal.  Even 
though  the  signal-to-noise  ratio  is  very  small,  and  thus  almost  lost  in 
noise,  the  phase  locked  loop  can  lock-in  on  only  a  few  cycles.  Laboratory 
experiments  with  the  565  PLL  chip  is  still  one  of  the  most  versatile  in 
many  applications  even  though  it  is  several  years  old.  It  has  a  bandwidth 
of  500  KHz,  a  center  frequency  stability  of  1200  PPM/°C  and  a  negligible 
frequency  shift  with  drift  of  the  supply  voltage.  A  circuit  diagram  for 
the  PLL  chip  used  is  shown  in  Fig.  C-2.  In  the  absence  of  signal,  the 
VCO  is  free  running  since  its  frequency  is  determined  by  external  compon¬ 
ents.  Fig.  C-3  shows  the  output  of  the  PLL  when  the  circuit  is  locked 
onto  the  received  signal  and  also  when  there  is  no  signal  present.  A 
filter  follows  the  PLL  in  order  to  eliminate  the  free  running  PLL  oscil¬ 
lator  signal.  The  particular  filter  used  in  the  laboratory  experiments 
is  shown  in  Fig.  C-4.  This  blocking  filter  is  followed  by  a  comparator 
so  as  to  return  to  an  appropriate  digital  form  of  the  signal.  An  active 
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filter  could  be  used  in  place  of  the  three  stage  passive  filter  shown 
in  Fig.  C-4. 

B.  Active  Filter: 

The  Active  Filter  (AF)  used  for  laboratory  comparison  with  the 
PLL  is  shown  in  Fig.  C-5.  The  gain  of  the  op-amps  was  10  and  the  Q 
was  10-.  The  laboratory  experiments  made  clear  that  a  PLL  is  preferred 
to  the  AF  for  the  following  reasons: 

1.  The  PLL  reached  lock-in  in  fewer  pulses  than  the  AF. 

2.  The  PLL  performed  over  a  greater  range  of  pulse  rate  drifting 
than  did  the  AF. 

Even  with  high  performance  ,op-amps  in  the  AF,  laboratory  results 
suggested  that  the  PLL  performs  more  effectively. 
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